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Abstract 

In this paper we study the Matovich-Pearson equations describing 
the process of glass fiber drawing. These equations may be viewed as a 
ID-reduction of the incompressible Navier-Stokes equations including 
free boundary, valid for the drawing of a long and thin glass fiber. We 



"This paper was partly written during the visit of Ph.L. and A.M. at Necas Center for 
Mathematical Modelling in March/ April 2008 and they wish to express their thanks 
for the hospitality during the stay. 

f The work of E.F. was supported by Grant 201 /09/0917 of GA CR as a part of the gen- 
eral research programme of the Academy of Sciences of the Czech Republic, Institutional 
Research Plan AV0Z10190503. 

"'"Research of A.M. on the mathematical modeling of glass fiber drawing was supported 
in part by the EZUS LYON 1 INGENIERIE during the period 2004 - 2007. 



1 



concentrate on the isothermal case without surface tension. Then the 
Matovich- Pearson equations represent a nonlinearly coupled system of 
an elliptic equation for the axial velocity and a hyperbolic transport 
equation for the fluid cross-sectional area. We first prove existence of 
a local solution, and, after constructing appropriate barrier functions, 
we deduce that the fluid radius is always strictly positive and that the 
local solution remains in the same regularity class. To the best of our 
knowledge, this is the first global existence and uniqueness result for 
this important system of equations. 

Keywords: Fiber drawing; Matovich-Pearson equations; incom- 
pressible free boundary Navier-Stokes equations; non-local transport 
equation; iterated comparison; barrier function. 

Mathematics Subject Classification MSC2000 (AMS class- 
code): 35Q30; 35Q35; 35R35; 35L80; 76D05; 76D27 

1 Introduction 

The drawing of continuous glass fibers is a widely used procedure. Industrial 
glass fibers are manufactured by a bushing with more than thousand noz- 
zles. Bushings are supplied with a molten glass from a melting furnace. Its 
temperature ranges from 1300-fT to 1800A'. In order to understand the glass 
fiber forming process, it is important to study the drawing of a single glass 
fiber. This is, of course, a significant simplification because we disregard 
interaction between fibers and between fibers and the surrounding air. For 
a single glass fiber, the hot glass melt is forced by gravity to flow through a 
die into air. After leaving the die, the molten glass forms a free liquid jet. 
It is cooled and attenuated as it proceeds through the air. Finally, the cold 
fiber is collected on a rotating drum. 

The molten glass can be considered as a Newtonian fluid and the process 
is described by the non-isothermal Navier-Stokes equations for a thermally 
dilatable but isochoric fluid. Since we deal with a free liquid jet, the problem 
is posed as a free boundary problem for the Navier-Stokes equations, coupled 
with the energy equation. We refer to [5] for detailed modeling and analysis 
of the equations describing the stationary flow inside the die. There are 
several models proposed to describe the various stages of the flow of a molten 
glass from the furnace to the winding spool: the slow flow in the die (the 
"first phase" of the drawing), the jet formation under rapid cooling (the 
"second phase"), and the terminal fiber profile (the "third phase") (see [6]). 

Since we consider long (their length is approximately 10 m) and thin 
(their diameter varies from 1 mm to 10 /xm) fibers, it is reasonable to apply 
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the lubrication approximation to the model equation. This approach yields 
good results, at least for flows far from the die exit and in the so-called 
"third phase" of the fiber drawing. 

A standard engineering model for the isothermal glass fiber drawing in 
the "third phase" is represented by the Matovich-Pearson equations. For an 
axially symmetric fiber with a straight central line, they read 

d t A + d x (vA) = 0; d x (3fi(T)A 3 x v) + d x (a(T)VX) = 0, (1) 

where A = A(t,x) is the cross-area of the fiber section, v = v(t,x) is 
the effective axial velocity, 3// is Trouton's viscosity, and a denotes the 
surface tension. As the coefficients \i and a depend on the temperature, it is 
necessary to take into account an equation for the temperature T = T(t, x). 

The original derivation of the system ([TJ is purely heuristic and obtained 
under the assumptions that: (i) the viscous forces dominate the inertial ones; 
(ii) the effect of the surface tension is balanced with the normal stress at the 
free boundary; (Hi) the heat conduction is small compared with the heat 
convection in the fiber; (iv) the fiber is almost straight, and all quantities 
are axially symmetric. We refer to the classical papers by Kase & Matsuo 
[161117] . and Matovich & Pearson [19] for more details concerning the model. 

Another derivation of the model based on a lubrication type asymptotic 
expansion can be found in the work by Schultz et al. [9j [23] , Dewynne et 
al. [3 S] , and Hagen [12] , with more emphasis on the mathematical aspects 
of the problem. The (formal) asymptotic expansion is developed with re- 
spect to a small parameter e, proportional to the ratio of the characteristic 
thickness Re in the radial direction and the characteristic axial length of the 
fiber L. The fact that the viscosity changes over several orders of magnitude 
is surprisingly ignored in these studies. As a matter of fact, the viscosity 
coefficient depends effectively on the temperature, with values varying from 
10 to 10 12 Pa sec, while in the above mentioned asymptotic expansions it 
is considered simply of order one. A correct formal derivation was given in 
[T], and it is in full agreement with the model announced in [12] . Finally, a 
full non-stationary model of a thermally dilatable molten glass, with density 
depending on the temperature, was derived in [6]. 

A mathematical analysis of generalized stationary Matovich-Pearson 
equations is performed in [1] (see also [2]). The non-stationary case, without 
surface tension and with advection equation for the temperature, is studied 
by Hagen & Renardy [10] . They prove a local-in-time existence result in 
the class of smooth solutions. Their approach is based on a precise analysis 
of the dependence of the solution of the mass conservation equation on the 
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velocity. This method requires controlling higher order Sobolev norms in 
the construction of solutions by means of an iterative procedure and works 
only for short time intervals. Hagen et al. \\.\\ [T3| [T4] have also undertaken 
a detailed study of the linearized equations of forced elongation. Despite 
this considerable effort, global- in-time solvability of the Matovich- Pearson 
equations was left open. 

The ratio a/fi is small, and, furthermore, the inertia and gravity effects 
are negligible in most applications. Accordingly, we consider the Matovich- 
Pearson equations ([I]) with a = 0, meaning, the isothermal drawing with 
constant positive viscosity and in the absence of surface tension. For a pre- 
scribed velocity at the fiber end points, the important parameter is the draw 
ratio, being the ratio between the outlet and inlet fluid velocities. It is well 
known that the instability known as a draw resonance occurs at draw ratios 
in excess of about 20.2. Linear stability analysis was rigorously undertaken 
by Renardy [22]. Moreover, in [26] . it was established that the cross section, 
given by the Matovich-Pearson equations with a = 0, may vary chaotically 
at a draw ratio higher than 30, under the condition of periodic variations 
of the input cross section. There are also numerous articles devoted to 
numerical simulations confirming such a conclusion. Fairly complete simu- 
lations can be found in the papers by Gregory Forest & Zhou [HI [25]. Their 
simulations predict various aspects of the physical process, like a linearized 
stability principle, bounds on the domain of convergence for linearly stable 
solutions, and transition to instability. Their analysis completes that of [7]. 

The above mentioned simplification of system ([I]) is briefly discussed in 
[15] , however, without any rigorous proofs. Our idea is to use the particu- 
lar structure of the system with a = 0, and to prove short-time existence 
of smooth solutions satisfying good uniform estimates. Then, performing 
a qualitative analysis of the solutions and constructing appropriate barrier 
functions, we show that the cross-section area remains bounded below away 
from zero. This observation allows us to deduce existence as well as unique- 
ness of global- in-time solutions. 

2 Isothermal fiber drawing without surface ten- 
sion 

We study the system of equations 





(3) 
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supplemented with the boundary and initial conditions 

A(t,0) = S o (t) in (0,T), A(0,x) = S 1 (x) in (0,L), (4) 
v(t,0)=v in {t) in (0,T), «(*,L)=«i(*) in (0,T). (5) 

Here v is the axial velocity and A denotes the cross section, L, T are given 
positive numbers, and 3/U > denotes Trouton's viscosity assumed to be 
constant. 

The data satisfy 

{0 < v m < v in (t) < v L (t) < V M for any t G (0,T), 
(6) 
< S m < 5 (t),5i(x) < S M for all (t,x) G Q T , 5 (0) = Si(0). 

Moreover, the functions Uj n , u^, 5o, 5 1 ! belong to certain regularity classes 
specified below. 

2.1 A priori bounds 

Our construction of global-in-time solutions is based on certain a priori 
estimates that hold, formally, for any smooth solution of problem ([2]) - ([5]), 
with the cross-section area A > 0. The crucial observation is that, as a 
direct consequence of (|3|), 

A{t,x) d x v(t,x) = X (t) for any t G (0,T), (7) 

where % is a function of the time variable only. Moreover, as A is positive 
and the axial velocity satisfies the boundary conditions ©, we deduce that 

x(t) > for any t G (0, T), (8) 

which in turn implies 

d x v(t,x) > for all (t,x) G Qt- (9) 

Accordingly, 

Vin(t) < v(t,x) < VL(t) for all (t,x) G Qt- 
Next, we rewrite equation ([2]) in the form 

d t A + vd x A = -x < 
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yielding 

A(t,x) < Sm for all (t,x) G Qt- (10) 

Integrating (J7J) over (0, L) and using ([5]) and (|10|) give rise to the uniform 
bound 

< X (t) < ^ ~ Vinit) S M for all t £ (0,T). (11) 

In order to deduce a lower bound for the cross section area A, we first 
observer that A and d x A satisfy the same transport equation, namely, 

d t A + d x {vA) =0, (12) 
d t (d x A) + d x (v(d x A)) = 0. (13) 

In particular, 

d t (d x \og(A)) + vd x (d x \og(A)) = 0. (14) 
Evaluating the boundary values of d x log(^4) with namely, 

(9 x log(^)(i,0) = - v , a x log(A)(0,x) - 



S (t) ' x m /v ' ' S 1 (x 
where, in accordance with ([2]), ©, ©, ©, and (HH) 

^ (t ' 0) = -sk ( x< ' ) + ^ (t) 

> ( ^C)-»-W SM + ^( t ) 

We deduce easily the desired lower bound on A 

A(t, x) > Am > for all (t, x) G Q T , (15) 

where the constant A m is determined solely in terms of v m , Vm, S m , Sm, 
and the first derivatives of So, S±. 

The a priori bounds derived in ((7|) - (|15p form a suitable platform for 
the existence theory developed in the remaining part of this paper. 

3 Short time existence of regularized strong solu- 
tions 

In addition to ([6]), we shall assume that 

5oG^ 2 '°°(0,T), S 1 eH 2 (0,L), v^vlEW^^T), (16) 
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where the symbol W k,p denotes the standard Sobolev space of functions 
having k— derivatives L p — integrable, and H 2 = W 2,2 . For further use, we 
introduce the notation 



Q°>° = -^(0)-v in (0)^(0). (17) 

Let us begin with a list of definitions: 

Definition 1. Let to be a positive number. A pair (A,v), defined on Qt = 
(0,io) x (0, L), is a strong solution of (H|)-([5]) if 

AeW l ^((0,t )x(0,L)), (18) 

v, 8 t v, d x v, d 2 x v, d 2 v E L°°((0,t ) x (0,L)), (19) 

(A,v) satisfy equations (fjj) - (0) a. e. in (0,t ) x (0,L), (20) 

A > on Q to and A satisfies Q pointwise , (21) 

v satisfies @ pointwise. (22) 

Definition 2. Let t be a constant, < t < T. For h G L°°(0, t ; H 1 ^, L)) 
with d x h E L°°(Q,L;L 2 (0,to)), we define the energy functional £ as 

£(h) 2 = ess sup \\h(t, Oll|i(o,i) + v m ess sup \\d x h(-,x)\\ 2 L 2 {0M) . (23) 

0<t<to 0<x<L 



Similarly, the radius R is defined by 

2' 



l R2 =L s Kf) \ {x)dx+ L \ sl+ k 



dS 2 



dt 



(t) dt. 
(24) 



Definition 3. For R given by \2J$ and < to <T, we denote by S(to,R) 
the convex set of nonnegative functions h defined on Qt such that 

h G W 1>oo (Q,t ;L 2 (0,L)) n L°°(0,t ; H 1 ^, L))n 

n^(0, L; L 2 (0, t )) n L°°(0, L; H\0, to)), (25) 

£(h) < R, and 

ess sup \\d t h(t,-)\\ L 2 ( on+ sup \\d t h(; x)\\ L 2, ,t ) < 
o<t<t xe[o,L] 



RV M \2 + \ + i^\ +\Q <°\(y/T + VZ) (26) 



L L 3 / 2 / 

h(0,x) = S 1 (x) and h(t,0) = S (t). (27) 
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Definition 4. For R given by {21$ , 

a = A| Q o >0|2 + ± Sm L\Q°'% f> = + & M V M + 

we denote t* £ (0, T] a positive time satisfying 

o r3/2 p2 

e< ^^_ and f< (28) 

In this section, for <5 > small enough, we construct a family of approxi- 
mate solutions (A s ,v s ) solving the following initial-boundary value problem: 

8 t A 5 + d x (v s A s ) = in Q 5 tQ = (5, t ) X (0, L), (29) 



Vl(S) - Vi n (5) t . 

— — - = m (0,5] x (0,L), (30) 



c^/M 5 c^ 5 ) = in Q t0 = (0, t ) x (0, L), (31) 

A 5 (t,0) = 5 (t) in (0,T), A a (0, a?) = Si (a?) in (0,L), (32) 

?/(t,0) = ^„(t) in (0,T), v 5 (t,L) = v L (t) in (0,T). (33) 

Specifically, we prove the following result. 



Theorem 1. Let VL,Vi n ,So and S± satisfy il6\) . Consider t* > given by 
H28\) , to £ (0,t*), and let 5 £ (0, min {1, to}) be a small number satisfying 

^-5\Q°>°\>0. (34) 



Then the initial-boundary value problem [29]) - [33]) possesses a unique 
solution (A s ,v s ) in the class 

A 5 G C\[0, t ]; ^(0, L)) n C([0, t ]; H 2 (0, L)), (35) 
v S e C\[0, t ];H 2 (0, L)) n C([0, t ]; # 3 (0, L)). (36) 

Note that for a strictly positive /i £ S(to,R), the corresponding velocity 
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field v solving d x (hd x v) = with boundary conditions ([5]) reads 



(37) 



Jo Kt,0 

Remark 1. Under the compatibility condition 




(38) 



at dx r L dx 




(39) 



we could set 5 = 0. However, imposing \39\) is not physically justified. Note 
that this condition is systematically used in 111 Of as well as 115^ - 

The rest of this section is devoted to the proof of Theorem [TJ The basic 
and rather standard idea is to construct a sequence approaching a fixed 
point of a suitable nonlinear mapping. The proof is carried over by means 
of several steps. We fix to £ (0,t*) and 5 S (0, min {1, to}) such that ([3l|) is 
satisfied. 



We take an arbitrary ^4° G S(to,R). Then we substitute h = A° into 
([37]) and calculate v = v. We note that A € 5(to,-R) implies 



d x v°(t,x) > and v in (t) < v°{t,x) < v L {t) in [0,t ] x [0,L]. (42) 
Next, we introduce functions Q° and Q 0,5 defined on [0,to] by 



STEP 1 



v° e W^QO, t ]; H\0, L)) n L°°([0, t ]; H 2 (0, L)), 



(40) 



o Ou n\ v L{t)-V in {t) 1 , 

t^d^foO) = — ^ — w -, , G (0,t ), 




(41) 




(43) 



Q°' s (t) 



= < 




(44) 
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Obviously, Q°> 5 E W 1 ' 00 ^,^), and, using Jensen's inequality and ©, we 
get 

0<Q°(t) < , VL ® < V -^[ L A%x) dx 



o 



(i 



A°(t,x) 



< |^ll^(V)lb M < («) 
since ^4° E S(to,R), as well as 

Q°^)|<IQ°'°I +^ (46) 

Now we are in a position to define the solution operator: 
For A E S(t ,R) and v° given by (|3ZD with A instead of /i, we solve 
the initial-boundary value problem 

d t u(t,x) = -v°(t,x)d x u(t,x) ~ Q° >5 (t), (t,x) E Q t0 , (47) 
ii(0, x) = Si(x), x E (0, L), u(t, 0) = So(t), t E (0, to)- (48) 



Because of (|40p -(|42p and regularity and compatibility of the data, we may 
apply a result of Hagen &: Renardy [10] (see Theorem [3] in Appendix) to 
problem (|47p -( j48j) to obtain a unique solution u E C 1 ([0, to]; H 1 ^, L)) n 
C([0,t ];i3" 2 (0,L)). We set 

A x (t,x) = u(t,x), (M)eQ v (49) 

Relation (|49p defines a nonlinear operator, assigning to a given A E 5 (to, -R) 
the unique function A 1 in the class C^QO, t ]; -ff^O, L))nC([0, t ];H 2 (0, L)). 
By (|4"T|) . dfu+dtQ°' S E L°°(Qt ). Using a compactness lemma of Aubin type 
(see e.g. |24j ) we conclude that our nonlinear operator, defined on 5(to, R), 
is compact. 

At this stage we follow the ideas of |10] and [15] to show that this non- 
linear operator has a fixed point for all to < t*. They study the nonisother- 
mal fiber spinning and their system of equations is different. Hence we are 
obliged to give an independent proof of short time existence, but the result 
remains close to their considerations. 

The natural approach is to apply Schauder's fixed point theorem. To use 
it we have to prove that S (to, R) is a relatively compact convex set invariant 
under our nonlinear mapping. Finally, we should establish the continuity of 
the mapping (|49|) . 
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STEP 2 



We establish uniform L°°— bounds for the function u. We have 

Lemma 1. Let R be given by Definition^ and let t$ < t* . Then the solution 
u of problem Mty-fiiJtfy satisfies the estimate 

q 29 

-f < u(t, x) < S M + -j 1 < 2S M in [0, t Q ] x [0, L]. (50) 
Proof. Let s m G W 2 '°°(0, to) be the solution to the Cauchy problem 

^ = -Q°> 5 (t) in (0,td), Sm (0) = (51) 
at 3 

Then s m clearly solves (|47p and it follows from ([6]) and ()48|) that 

u(0, x) = S'i(x) > 5 m > s m (0), x £ (0, L), 
while ([6]) , (j34"l) , the nonnegativity of Q° , (flHj) , and (l5Tj) ensure that 

9 C 9 C /-mill {<5,i} 

" JO "J Jo 

9C /■min{5,t} , ,> 90 X 

-3/0 Q°'°{ 1 -'i)^< -y 1 + \Q % 

< S m <S (t) = u(t,0). 

The comparison principle then entails that u(t,x) > s m (t) for (t,x) £ Q to , 
whence, bv flggD , (p4")l . and (Hp) , 

9 c 

> Sm (t)>^i - / |Q°' 5 ( s )|ds 



25 m _ i n o,Oi^ _ VmR ,* . 5V„ 
" 3 N l 2 L3/2 " 4 • 

This proves the lower bound. Proving the upper bound is analogous. The 
comparison function is now given by the solution to the Cauchy problem 

= -Q°> 5 (t) m (0,t ), SM (0) = S M + ^. (52) 

Indeed, sm clearly solves (|4T|l with sm(0) > Sm > S\(x) = u(0,x) for 
x e (0,L) and it follows from ©, ®, ®, (06]), gH}, and ([52]) that 

aw(t) > s M (0) ~ |Q°'°|| - ^ t* > SW > u(t,0). 
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Applying again the comparison principle gives u(t,x) < for (t,x) G 

Qt , which completes the proof since (fM|) guarantees that 

s M (t) <S M + ^ + \Q°>°\ 5 - < S M + 



This proves the Lemma. 



□ 



Now we use equation (|37|) to calculate G(t) = d x u(t,0) G C([0,to]) 
getting: 



G(t) = a a .«(t,o) 



i 



d5o(*) 



+ Q°' s (t)\. 



Vin{t) V ^ 

Due to the assumptions on the data, G G VF 1,oo (0, to) an d we have 



|G(*)I < 



1 



Vin(t) 



dS (t) 
dt 



\Q°> s (t)\ 



(53) 



(54) 



Next, we take the derivative of equation (|47p with respect to the x vari- 
able. This yields that S = d x u solves 



d t S = -v°d x S - d x v°S, (t, x) G Q t0 , 



(55) 



5(0, x) = ^(x), x G (0, L); S(t, 0) = G(t), t G (0, t ), (56) 
ax 

the function v° being still given by (J37J) with A instead of h. 

Multiplying equation (|4T|) by u, equation ([55]) by S, integrating both 
equations on (0, t) x (0, x), (t,x) G Qt , and adding the obtained identities, 
we deduce: 

\ £ {u 2 + S 2 } (i, Odt + l J* v°(r, x){u 2 + S 2 }(t, x) dr+ 



t rx 
JO 



t rx 



Q U '°(r)u(T,0 d&T 



d x v (T^){S 2 -u 2 }(r,C) didr + 

\[{ si+ (^) 2 } {i)di+ \ii vm{r) ^ + ^ (r) dT - (57) 



Now we use 



Sb with (I53D to get 



V ln (T)G 2 (T) dT < 
2 



Vin(r) 



dS { 



dr 



o v in (r) 
dr + 



dT 



dS ° +IQ ' 5 (r)| 2 | dT< 
4min{M} ^ 
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The third term on the left hand side of (j57[) is estimated with the help of 
H and (EDI) as 



f f d x v°(r,0(u 2 - S 2 )(t,0 didr < f f Q^Lu 2 (t,0 d^dr 
JO JO Jo Jo A ( T ;£) 



dT<AS 2 M / (v°(T,x)-v in (T)) dr 
Jo 



< 4S 2 M V M t, 



(59) 



while the fourth term obeys 
ft ex 

i0,<5 



JO 



Q u ' d (r)u(T,C) d£dT 



<2S M L f \Q°' S (t)\ dr 
Jo 



< 2S M L ( \Q°<°\ + 



VmR \ 

L 3 / 2 J 



t, 



(60) 



thanks to (|46D and (1501). 

Now, setting 

y(t)= sup { [ X {u 2 + S 2 )(t,0 d£ + f v\t,x){u 2 + S 2 }(t,x) dr) , 
xe[o,L] I Jo Jo J 

we may insert (f58l) - (f60l) into (jSTjl and use and (|28j) to get the following 
estimate: 



y(t) < 



/>+(t) 2 H 



R 2 R 2 t R 2 

< 1 < — . 

~ 8 8 t* ~ 4 



(t) dr + (a + fJ)t 



Recalling ([6]) and (J42J) , we thus conclude that 

8{u) < R. (61) 
Next, we use equation (j4"T|) and estimate (j6T1) to control dtu, obtaining 
sup \\d t u(t)\\ L 2 {0tL) + sup ||<9tu(-,x)|| L 2 (0>to) < 

t£[0,t ] x€[Q,L] 

RV M {2 + \ + ^) + |Q°'°|(V^ + VT). (62) 
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Therefore the image of a nonnegative function from S(to,R) remains in 
S(to,R) and satisfies the L°°— bound (|50p. Therefore our nonlinear map 
([4~9j) maps the convex set 

S (t , R) = {/ G S(t , R) | / satisfies ([50]) a.e.} (63) 

into itself. 
STEP 3" 

Let X be the intersection of the Banach spaces W 1,00 ([0,to]; L 2 (0, L)), 
L°°([0,t ];H l (0,L)), W^°°(0,L;L 2 (0,to)) and L°°(0, L; H l (0, t )). Clearly, 
5o(io 5 -R) is a convex, bounded and closed subset of the Banach space X. 
We apply the Schauder fixed point theorem to prove that the mapping (|49l) 
admits a fixed point in So(to, R). After Step 2, it remains only to prove the 
sequential continuity of the map (1491) . Hence let {A fc } fceN , be a sequence 
converging in So(to, R) to A. Then we have 



A k ^A strongly in C([0, to]; L (0, L)) n L (0, t ;C([0, L])), (64) 
Ql = l L ^Q°= Jl "J" UnifOTml y in C «°< *0]), (65) 



A k {t,x) J A(t,x) 

strongly in C([0, t ]; ^(0, L)). (66) 



Let u k be the solution to (|47|) - (|48|) . corresponding to A k and u the so- 
lution corresponding to A. Then by analogous calculations to those per- 
formed in Step 2, we get u k -»• it in C([0, t ]; fl^O, L)) (1L 2 (0, t ; C x (% L})). 
Using the equation (|47p. we find out that one also has <9tu fc — > dtu in 
C([0,t ];^ 2 (0,L)) nL 2 (0,t ;C([0,L])). Therefore the mapping flMD is con- 
tinuous and compact and by Schauder's fixed point theorem there is at least 
one fixed point A 5 G So{to,R)- 

Denoting the corresponding velocity field by v s , we have 

and since t h-> d^i.O) G M /1 ' oo (0, t ), we may apply Theorem [3] in Ap- 
pendix to conclude that, in fact, 

A 5 G CHIO, to]; ^(0, L)) n C([0, t ]; # 2 (0, L)). 
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It remains to prove uniqueness. 



STEP 4 



With the obtained smoothness, uniqueness is easy to establish. It suffices 
to notice that L 2 (0, to', H 1 ^, L)) perturbation of v is controlled by the L 2 - 
perturbation of A s in x and t. Then we use this observation, regularity of 
A s , and GronwalPs lemma to obtain uniqueness. 

The proof of Theorem [1] is now complete. 

4 Global existence of regularized strong solutions 

Now we suppose that the regularity of the solution or/and the strict posi- 
tivity of A s , stated in ()50|) . breaks down at the time t p . Our goal is to prove 
that t p = T. 

Theorem 2. Under the hypotheses of TheoremU\ the initial-boundary value 
problem i29\) - ^^) has a unique solution 



The remaining part of this section is devoted to the proof of Theorem [2j 



A S G C\[0, T];H\0, L)) n C([0, T];H 2 (0, L)), 
v S G C\[Q, T];H 2 (0, L)) n C([0, T];H 3 (0, L)). 



STEP 1 



First we recall that, owing to f)31 1) and (|33p . v 5 is given by 




(67) 




(68) 



and we put 




(69) 
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We point out that the upper bound in ([50]) is independent of t*. Indeed, 
we prove now that it is valid regardless the length of the time interval. To 
this end, we introduce the solution S 5 M G W 2 '°°(0, to) to the Cauchy problem 



dS s M (t) 
dt 



IQ°' |X[0,5](*) ; * 6(0, to), S° M (0) = S M . 



(70) 



Owing to the positivity ® of Q 5 , S 8 M is a supersolution to ([29]) - ([30]) with 
S 5 M {0) = 5 A / > Si(s) = A s (0,x) for x G [0, L] and S s M (t) > S S M (0) > 
S (t) = A 5 (t,0) for t G (0,t ). The comparison principle then implies 



A s (t,x)<S s M (t)<S M + S\Q 



0,0i 



(71) 



This proves the upper bound and, in addition, the estimate is independent 
of the length of the time interval. 

Therefore, by Jensen's inequality, we may infer that 
< Q 5 (t) < 



V M 



V M , 



'^r-(S M + 6\Q Q '°\), 



\G s (t)\ = \d x A\t, 



Vin(t) 



dS Q (t) 
dt 



+ Q cut > s (t) 



< 



dS (t) 



dt 



+ \Q 



0.0 



+ Yk {Sm + 5\Q^\) 



(72) 



(73) 



for any t < t p , where t p is the critical time at which A s (t, x) attains zero for 
some x. 



STEP 2 



Having established that d x A s (t, 0) is bounded in L°°(0,t) independently 
of t, we now look for a strictly positive lower bound for A s . 

Proposition 1. There are constants C\ and C2, independent of the length 
of the time interval and of 5, such that we have 



d< ft, log A 6 \t, x)<C 2 on Q to . 
Proof. We notice that y = d x log A s satisfies the equation 



d t y + v°d x y +js(Q°- Q cut > & ) = in Q tp , 



(74) 
(75) 
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and 
2/(0, x) 



d log Si, . . rN , , 

-(a;), xG(0,L), i/(t,0)- 



r/.r 



S (t) 



t£(0,t p ). (76) 



The function F(t,x) = (Q s - Q cut ' 5 )(t)/A s (t,x) vanishes for t > 5 and, due 
to ([50]) . satisfies the following bound 



|F(r) 



Il°°(o,l) 



min {t,<5} 



1 1^(^)11^(0^) dr 

0,0, , 



< ^-(IQ^I + IIQIIl-^) 



for every t since <5 < to- 

We next introduce the solution y m to the ordinary differential equation 



dy v 



dt 



■(t) = \\F(t)\\ LO o {0tL) y m (t) in (0,T), 



with initial condition 



^ mm | dl °^ x Sl } 



(78) 



(79) 



(80) 



Owing to ([73]) . ([76]) . ([To]) , and (j79j) . y m is nonpositive and thus a subsolution 
to ([75]) . and satisfies y m (i) < y m (0) < y(i, 0) for £ E (0, i p ) and y m (0) < 
y(0,x) for x 6 (0, L). The comparison principle then entails that 



with 



Gm = — 



dS 



dt 



L°°(0,T) 



+ | Q o,o | + ^ (5m + 5|q o,o | 



2/(M) > y m (t) = y m (0)exp / \\F(T)\\ Lao ^ L ) dr. 

Jo 

Since y m (0) < 0, we deduce from ([77]) that y(t,x) > C\ for some positive 
constant Ci, independent of (5 and io- 

Similarly, let Ym be the solution to the ordinary differential equation 



dYu 
dt 



(t) = \\F(t)\\ Loom Y M (t) in (0,T), 



(81) 
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with initial condition 

Y M (0) = max j max j ^jf 1 } , G M J > 0, (82) 

the constant Gm being defined in (f80|) . It follows from (f73|) . (f76|) . ([ST]) , and 
([82]) that Ym is a supersolution to (|75|) which satisfies Ym(£) > Ym(0) > 
y(t,0) for i £ (0,i p ) and Ym(0) > y(0,x) for x E (0, L). Using once more 
the comparison principle and (|77h . we conclude that 

y(t,x) < Y M (t) = y A /(0)exp / \\F(r)\\ L oo i0>L) dr < C 2 , 

Jo 

the constant C2 being independent on 5 and t. □ 
Since log^ 5 (t,0) = log So (i) > logS^ > -00 by © and 
log A* (t,x)= log S (t) + [ X d^ogA 5 (t,0 de, 

JO 

we infer from (|74p that 

II log A <5 || L oo( 0jto;H /i,oc,( 0jL )) < C. (83) 

Therefore is strictly positive on Q t for all t < t p , and bounded from 
below and from above by a constant which is independent of both t and 5. 



STEP 3 



Having established (|83p we easily obtain the following estimates: 

||^IU~( Qt) + ||^ 5 ||^ (Qt) <C, (84) 

\\v 5 \\L°°(Q t ) + \\dxV S \\ L °°(Qt) + \\dlv s \\ L oo {Qt) + ||^t«*|Uoo ( Q t ) < C, (85) 

where C is again independent of t and 5. 

The estimate (|85|) guarantees that the coefficients in equations ([29|) - ([30]) 
remain regular, whence Theorem 2] is applicable. Consequently, A s remains 
bounded in C^flO, t]; H 1 ^, L)) n C([0, t]; F 2 (0, L)). 

We conclude that, for all i, A s is bounded from below by a positive 
constant, independent of t, and the norm of A s in C 1 ([0, i]; -ff 1 (0, L)) n 
C([0, i]; H 2 (0, L)) remains bounded by a constant, also independent of t, 
that may, however, depend on 5. The maximal solution therefore extends 
to [0, T], and, in fact, we have established the existence of a unique strictly 
positive solution A s on (0, T) x (0, L). The corresponding velocity v s is given 
by (|37]), with h = A 8 . 

This completes the proof of Theorem [2] . 
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5 Existence of a unique strong solution 

At this stage, we are ready to establish the main result of the paper. 

Theorem 3. Under the hypotheses {TjJ), Mb}) , the initial-boundary value 
problem (H|)-([5|) possesses a unique (strong) solution A, v on Qt, belonging 
to the class 

A, d t A, d x AeL°°((0,L) x (0,T)), 
v, dtv, d£ x v, dl x v e L°°((0,L) x (0,T)). 

Proof. We just recall the estimates obtained in the proof of Theorem [21 valid 
independently of 5: 



< Ci < A\t, x)<C in Q T , (86) 
\l°°(Qt) + \\dtA 5 \\ L c (Q T ~ ) 



d x A s \\ Loo{QT) + \\dtA s \\ Loo{QT) < C. (87) 



Since L°°{Qt) is the dual space of the separable Banach space L 1 (Qt), 
Alaoglu's weak* compactness theorem gives weak* sequential compactness. 
Therefore there exist A and v such that 

A 5 ->■ A uniformly in Q^, as 8 ->■ 0, (88) 

a^^d^ weakly* in L°°(Q T ), as (5^0, (89) 

3^ d t A weakly* in L°°(Q T ), as 5^0, (90) 

5 ^ m , VL(t)-Vin(t) T ^ ., , - 7=r x _^ n 
« -> v = Uin(tJ H -T / uniformly m Q T , as d ->• 0. 

(91) 

d^" 5 -> 8 x v = VL ® ZjW 1 uniformly in Q^, as 5^ 0, (92) 
rL A(t,x) 

Jo A(t,£) 

at least for suitable subsequences. Obviously, (A, v) solves the system ([2])- 
([5]). Moreover, by virtue of the interior regularity, the equations (J2J)- ([3j) are 
satisfied pointwise. Finally, according to the smoothness of A and v, the 
proof of uniqueness is straightforward. □ 



Remark 2. The standard way of proving uniqueness relies on Gronwall's 
inequality. For any bounded time interval, small L? -perturbations of the 
data vl, %, So in the L? — norm result in the corresponding variation 
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of the solution in the same norm, that may depend exponentially on the 
length of the time interval. Better estimates would require refined analytical 
arguments. 

6 Appendix 

Here we recall the result from [10] , which is used in this paper: 

Theorem 4. Let f and p be given continuous functions defined on [0, to] x 
[0,6], 6 > 0. Let u° G C([0, b]) and u b G C([0,t ]) and let us supposqj that 

p,f€ W l >°°([0, t o ];H\0, b)) n L°°([0, t }; H 2 (0, b)), 

are such that d x p(-, 0) and d x f(-, 0) G i? 1 (0, to); (93) 

p < on [0,t ] x [0,6]; (94) 

u°G# 2 (0,6), u b € H 2 {0,t ); (95) 

u°(0) = u fe (0); «V(0) = p(0,0)a s u°(0) + /(0,0). (96) 

T/ien t/je boundary-initial value problem 

d t u = p(t, x)d x u + f(t, x), {x, t) G (0, t ) x (0, 6), (97) 
u(0, x) = u°(x), x G (0, 6); u(t, 0) = u b (t), t G (0, to)- (98) 

/ios a solution 

nG C 1 ([0,to]; J ff 1 (0,6))nC([0,t ]; J H' 2 (0,6)), (99) 
w/wc/i is untgue in W 1 ' 00 ^, t ; L 2 (0, 6)) n L°°(0, t ; F^O, 6)). 
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